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The quantum Ising model: finite sums and hyperbolic functions 
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Jagiellonian University, Institute of Physics, Lojasiewicza 11, 30-348 Krakow, Poland 

We derive exact closed-form expressions for several sums leading to hyperbolic functions and dis¬ 
cuss their applicability for studies of finite-size Ising spin chains. We show how they immediately lead 
to closed-form expressions for both fidelity susceptibility characterizing the quantum critical point 
and the coefficients of the counterdiabatic Hamiltonian enabling arbitrarily quick adiabatic driving 
of the system. Our results generalize and extend the sums presented in the popular Gradshteyn and 
Ryzhik Table of Integrals, Series, and Products. 


Introduction 


In typical condensed matter setups the size of a system is thermodynamic and so the finite-size effects are n^hgible. 


The situation can be quite the opposite with quantum simulators, i.e., quantum systems such as cold atoms [l|, and 
ions ii that are tailored to emulate the challenging condensed matter models. One of the ultimate goals of this 
field is to emulate spin systems that are analytically and numerically unapproachable in the thermodynamic limit. 
Realization of such a goal requires substantial experimental progress, but it is reasonable to expect that on-demand 
emulation of arbitrary spin models is within experimental reach | 
manage to emulate spin systems composed of less than 20 spins 

It is thus desirable to develop better understanding of finite-size effects in many-body spin systems to be able 
to (i) extrapolate the finite-size data to the thermodynamic limit; (ii) understand the limitations of the current 
generation of quantum emulators. Actually, the finite-size effects are interesting on their own in systems undergoing 
a quantum phase transition. Indeed, such systems are especially sensitive to the system size because their correlation 
length diverges around the critical point in the thermodynamic limit 00. Additionally, we mention that finite-size 
systems are relevant for quantum engineering, where the goal is to prepare a system in a desired quantum state (see 
e.g. Ref. (Tl|l. 

Due to its exact solvability, the one dimensional quantum Ising model is a perfect system for analytical studies of a 
quantum phase transition |9l.ll0l|. Such studies have been extensively performed in the past in the thermodynamic limit. 
This limit allowed for the observation of singularities at the critical point and simplified the analytical computations 


LTlie current state-of-the-art experiments, however, 
where the finite-size effects come into play. 


The calculations in the thermodynamic limit of the quantum Ising model rely on the replacement of sums with 
integrals. If we would be able to find closed-form expressions for such sums, we would get detailed analytical insights 
into the finite-size effects in a model system undergoing a quantum phase transition. The main goal of this work 
is to provide several closed-form expressions for the sums that appear in the studies of the quantum Ising model. 
Additionally, we briefly discuss two examples, where the use of our sums immediately leads to elegant expressions 
for both fidelity susceptibility capturing the key properties of the quantum critical point and the coefficients of the 
counterdiabatic Hamiltonian allowing for implementation of perfectly adiabatic dynamics on an arbitrarily short time 
scale. 


Results 


There are four interestingly looking sums in Chapter 1.382 of Ref. 
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which are valid for integer n > 1 (we assume such n everywhere in this article). To explain where these sums come 
from, we consider polynomials 
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which can be also written as 
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If we now substitute 


g = exp(a;) 


(9) 


into Eqs. (H])-® and differentiate them with respect to x, we will get Eqs. after simple algebra. As will be 

explained in the Discussion section, the parameter 5 ® is the magnetic field in the quantum Ising model. 

The sums o-@ can be generalized in the following way 
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Eqs. (I10|) - (|lll) are valid for m = 0,1, • • • , 2n, while Eqs. (P^-([TH|) hold for m = 0,1,-- - ,2n+ 1. These results 
reproduce Eqs. for m = 0. 

Additionally, we propose another set of identities 
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Eqs. (fTH) and (ITS)) yield a correct result for m = 1, 2, • • • , 2n — 1, while Eqs. (fTCl) -(fT7 |) hold for to = 1, 2, • • • , 2n. 

The sums (nni)-® have been conjectured by applying the techniques described in Appendix A of Ref. [l^ . Two of 
them, (nni) and m, have been already proven in Ref. |lq . We will outline below the steps allowing for an inductive 
proof of all the sums. This way we will generalize the former results and simplify a bit the proof. 

To compactify notation, we define two functions 
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where the momenta k are chosen such that replaces in Eqs. (IT0l) - (|17p . Next, we show that 
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where A^ = ‘2’^^ cos{mk). These recurrence relations can be derived by substituting 
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and 


sin^(fc) =2[cosh(a;) + cos(fc)][sin^(fc/2) + sinh^(a:/2)] — sinh^(a;) 


( 22 ) 


into cos[(to + l)fc] = cos{mk) cos(/c)—sin(mfc) sin(fc) and sin(fc) sin[(TO + l)fc] = sin(fc) sin(TOfc) cos(fc)+cos(mfc) sin^(fc). 
To proceed with the inductive proof, we find in an elementary way that 
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for m = 1, 2, • • • , 2n — 1 and 
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for m = 1, 2, • • • , 2n. 

First, we assume that fm{x) and hm{x) are given by the right-hand sides of Eqs. (Unil-dlSl) and (Hi-dni) , respec¬ 
tively. 

Second, we check that our inductive conjecture holds for m = 1. This is done by substituting Eqs. m and (1^^ 
into Eqs. m and (ITO)) taken at m = 1 and then comparing the resulting expressions to Eqs. (HD-®. 

Third, it is a straightforward exercise to check that the expressions for fm+iix) and hm+i{x) obtained from the 
recurrence equations (l20ll are exactly the same as the right-hand sides of the conjectured equations (fTOll -dTT i after 
the replacement of m by m + 1. The last remark is true provided that expressions for A^ are used in the recurrence 
equations for appropriate m listed below Eqs. (l23l) and (l24l) . This completes the proof that Eqs. (uni), (HI]), dni), and 
(fT^ are correct for m = 1, 2, • • • , 2n — 1, while Eqs. m, uni, dnj), and (flTl) hold for m = 1,2, • • • , 2n. Additionally, 
one may substitute m = 2n into Eqs. uni) and m and m = 2n + 1 into Eqs. m and m to immediately verify 
with Eqs. 0-® that Eqs. (nnD-Gi are correct for such m as well. This concludes the proof of Eqs. (Iin|)-(1I7|). 


Discussion 

The results obtained in the previous section provide useful analytical tools for the studies of the quantum Ising 
model. The Hamiltonian of this model reads 


N 

^ = -E(«+i+5^*)’ (25) 

i=l 

where N is the number of spins, ai denote the Pauli matrices acting on the i-th spin, g is a magnetic field, and 
periodic boundary conditions are assumed = dj- This exactly solvable model has many interesting properties 

SE3. We will mention two of them that are relevant for the present discussion. 

First, the above Hamiltonian commutes with the so-called parity operator 

N 

p=n-.' 

i=l 

Therefore, its Hilbert space can be divided into positive and negative parity subspaces (eigenvalues of the parity 
operator are ±1 because = 1). The parity of the ground state depends on both the system size and the direction 
of the magnetic field (see e.g. Ref. [13). 

Second, the quantum Ising model can be mapped to the non-interacting lattice fermionic model through the Jordan- 
Wigner transformation and then it can be diagonalized in the momentum space [l^, [T^ . There are four possibilities 
for the choice of momenta fc in a periodic system (see e.g. Ref. ca). In the positive parity subspace of the even-sized 
system 
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while in the negative parity subspace of such a system 


A: = 0, ± —,7r, j = 1, 2, • • • , n — 1, N = 2n. 


In the positive parity subspace of the odd-sized system 


2n -|- 1 

while in the negative parity subspace of such a system 
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Forgetting about the irrelevant sign of k and the 0 and tt momentum modes, the summations over these four choices 
of momenta correspond to the four summations over j seen in Eqs. (|10l) -(fT3 l) and (fT^ -lfT ^ . 

What’s more, the denominator of all the sums that we have computed, sin^(fc/2) -|- sinh^(a;/2), can be mapped to 
the form familiar from the Ising model computations, 

— 2g cos{k) + 1, (30) 


after setting 


X = ln(5). 


(31) 


Note that we have already used this mapping at the beginning of the Results section Since both the eigenvalues 
and eigenstates of quantum Ising model depend on the factor (1501) - see e.g. Ref. Q - it is not surprising that the 
sums involving it appear in the studies of this model. 

Finally, we note that the sums involving positive integer powers of the factor (1301) in the denominator can be 
easily generated through the differentiation of the sums (iini)-(ini) with respect to x. We are now ready to show two 
applications for the results of the previous section. 

Fidelity approach to quantum phase transitions. The location of the critical point and the critical exponent 
quantifying power-law divergence of the correlation length around the critical point can be extracted from the overlap 
between ground states corresponding to different external fields [l9l - l^ . Such an overlap is called fidelity [a. It 
defines the so-called fidelity susceptibility y through a Taylor expansion in the field shift S [1^ 

\{g\g + S)\ = l-xi 9 )^-^ + 0{6^), (32) 

where \g) is a ground state of a Hamiltonian depending on an external field g. It can be easily shown that [l^ 


xig) 


1 V ^ sin (fc) 
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in the quantum Ising model ((25l) : the restriction to 0 < fc < tt in the sum is applied to momenta from Eqs. (I26l) -(l29 l) . 

Therefore, to compute the exact closed-form expression for fidelity susceptibility xig) one needs to (i) determine the 
parity of the ground state and choose the appropriate momenta k from the list (I26l) -(l29 l) : (ii) choose the relevant sum 
from Eqs. dH-dUl), set TO = 1 to obtain sin^(A:) in the numerator of the summand, and compute derivative of the 
resulting expression with respect to x to square the summand’s denominator; (iii) substitute (I3I|) into the resulting 
expression. One gets in the end 


1652(|g|jv 2)^ 16g^{\g\^ + l){g^ - 1) 


(34) 


for all system sizes and magnetic fields. This expression was computed in a bit more complicated way in Refs. (T7l.[^. 

The quantum Ising model has the critical point at gc = 1 separating the ferromagnetic phase (0 < g < I) from 
the paramagnetic phase ((/>!). It is thus expected that its ground states are most sensitive to the change of the 
magnetic field near the critical point. Therefore, fidelity (1321) should have minimum near the critical point, and so 
fidelity susceptibility should have maximum there. One can use Eq. p4l) to easily analytically compute the location 
of the maximum of fidelity susceptibility and quantify how fast it approaches gc = 1 as the size of the system is 
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increased. One can also study with the help of Eq. ((Ml) the dependence of fidelity susceptibility on the system size 
and the distance from the critical point, which yields further insights into the quantum criticality of the Ising model 
(see Ref. jl^ for the details). 

Counterdiabatic driving. Perfectly adiabatic evolutions of a quantum system subjected to arbitrary time vari¬ 
ation of an external field can be engineered by a proper modification of a Hamiltonian (ill. |27H^ . 

For the quantum Ising chain this can be obtained after adding to the Hamiltonian (l25)) 
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describes interactions of m-1-1 adjacent spins. This so-called counterdiabatic Hamiltonian was computed for even-sized 
chains evolving in the positive parity subspace in Ref. [s^. We have checked that it can be written in such a form 
for systems of an arbitrary size evolving in a subspace of arbitrary parity. 

It should be perhaps explained that the fairly complicated counterdiabatic Hamiltonian is obtained by “reverse¬ 
engineering”. It is written down first in the momentum space, where it takes a simple form, and then it is transformed 
to the real space form. It is thus assumed that the system evolves in the subspace of the Hilbert space of definite 
parity. Such an assumption can be made because parity is conserved during time evolutions with the counterdiabatic 
Hamiltonian: [+ -f AH, P] = 0. 

The strength of the (m -I- l)-body interactions in the counterdiabatic Hamiltonian is given by 


, . 1 sin(fc) sin(mfc) 

= ^ 2^ g2 _2gcosik) + l 

0<k<-K ^ ^ 


(35) 


Using mapping ([5]) as well as Eqs. (I26I) - (I29I) and (m-mi), one immediately finds that for any system size TV > 1 
and TO = 1, 2, • • • , — 1 

Sg^+Pl±g^)’ 

where the -I- (—) sign should be used for evolutions taking place in the positive (negative) parity subspaces of the 
Hilbert space. The closed-form expression for Eq. (I35|) in even-sized systems evolving in the positive parity subspace 
was previously reported in Ref. . 

Summarizing, we have derived several expressions for the sums leading to hyperbolic functions. They generalize 
and extend the results listed in the popular Gradshteyn and Ryzhik book [l^. We expect that they will be useful 
in different studies and illustrated their applicability to the quantum Ising model, which has provided the primary 
motivation for this work. 
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[1] Lewenstein, M. et. al. Ultracold atomic gases in optical lattices: mimicking condensed matter physics and beyond. Adv. 
Phys. 56, 243 (2007). 

[2] Simon, J. et. al. Quantum simulation of antiferromagnetic spin chains in an optical lattice. Nature 472, 307 (2011). 

[3] Bloch, L, Dalibard, J. & Nascimbene, S. Quantum simulations with ultracold quantum gases. Nature Phys. 8, 267 (2012). 

[4] Blatt, R. & Roos, C. F. Quantum simulations with trapped ions. Nature Phys. 8, 277 (2012). 





7 


[5] Korenblit, S. et al. Quantum simulation of spin models on an arbitrary lattice with trapped ions. New J. Phys. 14 , 095024 

( 2012 ). 

[6] Lanyon, B. P. et al. Universal digital quantum simulation with trapped ions. Science 334, 57 (2011). 

[7] Islam, R. et al. Emergence and frustration of magnetism with variable-range interactions in a quantum simulator. Science 
340 , 583 (2013). 

[8] Richerme, P. et. al. Experimental performance of a quantum simulator: Optimizing adiabatic evolution and identifying 
many-body ground states. Phys. Rev. A 88, 012334 (2013). 

[9] Sachdev, S. Quantum Phase Transitions (Cambridge University Press, 2011). 

[10] Sachdev, S. & and Keimer, B. Quantum criticality. Phys. Today 64, 29 (2011). 

[11] Torrontegui, E. et al. Shortcuts to Adiabaticity. Adv. At. Mol. Opt. Phys. 62 , 117 (2013). 

[12] Lieb, E., Schultz, T. & Mattis, D. Two soluble models of an antiferromagnetic chain. Ann. Phys. (N.Y.) 16 , 407 (1961). 

[13] Barouch, E. & McCoy, B. M. Statistical Mechanics of the XY Model. II. Spin-Correlation Functions. Phys. Rev. A 3 , 786 
(1971). 

[14] Pfeuty, P. The one-dimensional Ising model with a transverse field. Ann. Phys. 57 , 79 (1970). 

[15] Gradshteyn, I. S. & Ryzhik, I. M. Table of Integrals, Series, and Products, 7th ed. (Academic Press, San Diego, 2007). 

[16] Damski, B. Counterdiabatic driving of the quantum Ising model. J. Stat. Mech. (2014) P12019. 

[17] Damski, B. & Rams, M. M. Exact results for fidelity susceptibility of the quantum Ising model: the interplay between 
parity, system size, and magnetic field. J. Phys. A 47 , 025303 (2014). 

[18] Barouch, E., McCoy, B. M. & Dresden, M. Statistical Mechanics of the XY Model. I. Phys. Rev. A 2, 1075 (1970). 

[19] Zanardi, P. & Paunkovic, N. Ground state overlap and quantum phase transitions. Phys. Rev. E 74 , 031123 (2006). 

[20] Albuquerque, A. F., Alet, F., Sire, C. & Capponi, S. Quantum critical scaling of fidelity susceptibility. Phys. Rev. B 81 , 
064418 (2010). 

[21] Gritsev, V. and Polkovnikov, A. in Understanding in Quantum Phase Transitions edited by Carr, L. (Taylor & Francis, 
Boca Raton, 2010); e-print arXiv:0910.3692 (2009). 

[22] Zhou, H.-Q., Zhao, J.-H. & Li, B. Fidelity approach to quantum phase transitions: finite-size scaling for the quantum Ising 
model in a transverse field. J. Phys. A 41, 492002 (2008). 

[23] Rams, M. M. & Damski, B. Quantum fidelity in the thermodynamic limit. Phys. Rev. Lett. 106 , 055701 (2011). 

[24] Gu, S.-J. Fidelity approach to quantum phase transitions. Int. J. Mod. Phys. B 24 , 4371 (2010). 

[25] You, W.-L., Li, Y.-W. & Gu, S.-J. Fidelity, dynamic structure factor, and susceptibility in critical phenomena. Phys. Rev. 
E76, 022101 (2007). 

[26] Damski, B. Fidelity susceptibility of the quantum Ising model in a transverse field: The exact solution. Phys. Rev. E 87 , 
052131 (2013). 

[27] Demirplak, M. & Rice, S. A. Adiabatic population transfer with control fields. J. Phys. Chem. A 107 , 9937 (2003). 

[28] Berry, M. V. Transitionless quantum driving. J. Phys. A 42 , 365303 (2009). 

[29] Deffner, S., Jarzynski, C. & del Campo, A. Classical and quantum shortcuts to adiabaticity for scale-invariant driving. 
Phys. Rev. X A, 021013 (2014). 

[30] del Campo, A., Rams, M. M. & Zurek, W. H. Assisted finite-rate adiabatic passage across a quantum critical point: Exact 
solution for the quantum Ising model. Phys. Rev. Lett. 109 , 115703 (2012). 


